An Upper Bound on the Zero-Crossing 
Distribution* 

By NICHOLAS A. STRAKHOV and LUDWIK KURZf 

Let Q(T) equal the probability that a random process, x{t), does not 
cross the zero axis in a given interval of length T, A family of upper bounds 
for Q(T) is derived with only weak restrictions imposed on x(t) and it is 
shown that for gaussian random processes only one member of the family 
provides useful formulae. Specific results are obtained for x(t) representing 
a number of interesting random pi-ocesses. 

I. INTRODUCTION 

Let Q(T) equal the probability that a random process, x(t), does 
not cross the zero axis in a given interval of length T. The problem 
of determining Q(T) (and related functions) has important appli- 
cations in communications theory and has been investigated by 
many authors. 1 -" Reference 5 gives an extensive bibliography of 
most of the related work on this subject prior to 1962. Despite all 
this effort, Q(T) is known only when x(t) is a simple nongaussian 
process (such as a process whose zero-crossings obey the Poisson 
distribution) or a stationary gaussian zero-mean process with one 
of four explicit correlation functions. 5 - 6 Most of the rest of the results 
obtained are either approximate or form upper or lower bounds. 5 

In this paper, we develop a whole family of upper bounds on 
Q(T). For computational purposes, however, only one member of 
the family has been found to provide useful results for most cases 
of interest. 

II. DERIVATION OF AN UPPER BOUND ON q(t) 

Consider the transformation 

1 f T 

Zt - y J sgn [x(t)] dt (1) 



* An abbreviated version of this paper was presented at the Fifth Allerton 
Conference on System and Circuit Theory. Montieello, 111., October 4, 1967. 
t New York University. 
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where x(t) is a sample function of a stochastic process,* T is a fixed 
observation interval, and z T is a random variable defined by the 
stochastic integral (1). The function sgn[x(t)] is defined as 

f+1, x > 
sgn [x] = i 0, x = 
-1, x < 0. 
Since z T is a random variable, it has a cumulative distribution func- 
tion, P(z T ), associated with it. From (1), two properties of P(z T ) are 
immediately apparent, regardless of the statistics governing x(t): 

(i) P(z r ) =0 for z r < -1 

and (2) 

P(z T ) = 1 for z T > 1 

(w) lim [P(l + «) - P(l - e)] = QrXT) (3) 



and 

where 
and 



lim [P(-l + e) - P(-l - c)] = <?l(T) (4) 

t—0 

Q n {T) = Prob {*(<) i0 for 0^<ir| (5) 



Qj(T) - Prob {x(0 gO for 0g^T|. (6) 

Obviously, Q(T) as defined previously is related to the last two 
quantities by 

Q{T) = Qu(T) + Q L {T). (7) 

If x (t) is a symmetricf process, then 

Qu(T) = Q L {T) = \Q{T). (8) 

As a consequence of properties (i) and (it), P{z T ) can be rep- 
resented by 

P{z T ) = G(z T ) + Ql(T)u(z t + 1) + Qtr(T)tt(2 T - 1) (9) 

♦Throughout this paper, we assume that almost all sample functions of the 
stochastic pocess are continuous. Thus, (1), (5), and (6) are well defined. 

t The stochastic process x(t ) will be called symmetric if the probability meas- 
ures that govern it also govern the process — x(t). 
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where G(z T ) is continuous at z r = ±1 and 

[0 for x < 



u{x) = - 

[l for x ^ 0. 

We assume throughout that the function G(z T ) is not identically 
equal to zero. If it were, then it would be easy to show that Q(T) 
is known exactly; that is, Q(T) = 1. 

Next, consider the even-order moments of P(z T ), denoted by the 
Stieltjes integral 

g«l= f Zr dP(z r ) k = 0, 1, 2, • • • . (10) 

By substituting (9) into (10) one obtains 

q 2k = f z 2 T k dG(z T ) + Q (T) + Q L (T) k = 0, 1, 2, • • • . (11) 
J -i 

Neglecting the first term in the right side of (11) (which is always 
positive) and taking (7) into account leads to a family of upper 
bounds for Q (T) expressed by 

Q(T) ^ q 2k k = 0, 1, 2, • • • . (12) 

For k = 0, (12) reduces to the obvious result 

Q(T) ^ l. 

Before discussing the usefuless of the inequality (12), an expression 
for the moments will be derived. 
From its definition, (10) , q- lk can be expressed as 

q ik = E{z 2 1 k \ 

where E{-} denotes the expected value of the quantity enclosed in 
braces. Substitution for 2* from (1) results in 

f } • • - j o v(.h)y(Q ■ ■ ■ yihu) dt, dt, ■ ■ ■ dt* k \ (13) 

where 

y(U) = sgn \x(t t )] i = 1, 2, ■ • • , 2k. 
Interchanging the order of integration and expectation yields 

q 2 k = j,2i J J '" I ^*> ' l - ' " ' ' * 2 *) dtl ^ 2 " " ' dl 2k (14) 
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where 

R(l t , k , ■ ■ ■ , Uh) = E{y(t x )y(t2) ■ • ■ J/W). (15) 

We now make some remarks concerning the ordering of the family 
of inequalities expressed by (12) . 

Denote the first term in the right side of (11) by s 2t , or 

s„ - f z% k dG{z T ). (16) 

We next establish that s, k > s 2l , +2 (k = 0, 1, 2, . . .) which in turn 
establishes 

1 > q 2 > q< ■ ■ ■ > Q(T). (17) 

The former inequality follows directly from 



= J' zV dG(z T ) 
^ f z- k dG(z T ) 



— s 2 k 
with equality if, and only if, G{z T ) is of the form 

G(z r ) = Au(z r + 1) + Bu{z r - 1). (IS) 

Since G(z r ) is continuous at z T = ±1, equality is not possible and there- 
fore 

S 2 *+a < s-, k 
which, together with (11), establishes (17). We next establish the 
readily proven fact that 

lim s 2 , < e, e > (19) 

k-* oo 

and therefore, 

\imq 2 k = Q(T). 

A— .oo 

We begin by choosing an a(k ) > such that 

f ' 0> z% k ' dG(z T ) <|, e> 0. 

This can always be done because G{z T ) is continuous at z T = —1. 
Using the definition of s^,(16), obvious symmetry properties, and 
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the fact that 

/ z J r dG{z T ) < / 4** rfG(z r ) 

» - i J -i 

for each k > k , it follows immediately that 

/•l /.l-n(JU) 

lim / g 2 £ (lG{z T ) < e + lim / 2" dG(*r). 

j-_cd .'-I <-_ .a. •'-H-a(*o) 

Since the sequence of functions {«"}, A" = 0, 1, 2, ••• is uniformly 
convergent to zero on the interval [— 1 + a(k ), 1 — a(k )], the limit 
and the integral may be interchanged yielding (19). 

In light of (17) and (19), it appears that (12) should be evaluated 
for as large a value of A- as possible. For the special case when x(t) 
is a stationary gaussian random process (assumed to be zero-mean 
without loss of generality), it does not seem to be possible to evaluate 
q 2 k for k > 1 as evidenced by the following discussion. 

As shown by McFadden, 7 the quantity R(t x , t 2 , • • • , /„), defined 
in (15), is equal to the sum of some simple terms plus a quantity P„(r), 
which is defined as 

PJi) = (27r)-" /2 |r|-* f dx t ■■■ T dx n exp \ -\ £ r7/**,l 

■'0 ''O & i,i J 

where r is a covariance matrix with elements 

/•„• = r(U - t t ) = E\x{t,)x(t,)\, i, j = 1, • • • , n 
J r I is the determinant of r 

/! rJfXtXj is the quadratic form associated with the inverse of r 

and 

Xi = x(t ( ), i = 1, 2, • • • , n 

In other words, P n (r) is the probability that the n jointly distributed 
gaussian random variables, x(t t ) (i = 1, • • • , n) are all positive. 

As discussed by McFadden, 7 and even more thoroughly by Slepian, 5 
expressions for P n (r) have not been obtained in terms of elementary 
functions for n > 3. Because of this fact, it seems unlikely that an 
expression for (14) with k > 1 can be obtained for a general gaussian 
process, x(t). It should be pointed out, however, that an expression 
for q A has been derived 8 for p(r) = exp (—\t\), but without first evaluat- 
ing PJi). This result is not included because, for this correlation 
function, Q(T) is known exactly. 5 
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ni. APPLICATION TO GAUSSIAN RANDOM PROCESS 

Assume that x{t) in (1) is a stationary, zero-mean, gaussian random 
process, normalized so that P (0) = 1 where p(r) = E{x(t)x(t+r)}. 
The relationship (14) will be evaluated for the case k — 1, that is, 
for 

q2 = Y 2 f f R{U ' k) dU dts (20) 

where 

R{t, ,t 2 ) = E {sgn x{U) sgn x(t 2 ) \ . (21) 

The latter expression has been evaluated by many authors (see page 
58 of Lawson and Uhlenbeck's book, 9 for example) and the result is 

R(t lf t>) =-sin- 1 [ P a - t 2 )]. (22) 

IT 

Substituting (22) into (20) and making the obvious simplifications 
in integration results in 

q 2 =- [ (1 - u) sin -1 [p(Tu)] du, 

or, in light of (12), 

Q(T) ^ - f (1 - u) sin -1 [ P {Tu)] du. (23) 

IT J 

This result has been obtained by Slepian, 5 who states it as Theo- 
rem 5. # Slepian's proof, however, is long and complicated, as opposed 
to the simplicity of the proof given here. Furthermore, extensions to 
other cases can be obtained using the new method. 

IV. APPLICATION TO SINE WAVE PLUS GAUSSIAN RANDOM PROCESS 

We now turn to applying (12) with ft; = 1 to the case where 

xit) = w(t) + A cos (2tt# + v) (24) 

where 

w(t) is a stationary, zero-mean, gaussian random process with 
normalized correlation function, p(r), 
<p is a random phase constant uniformly distributed on [0, 2v] , 
/ is the sine-wave frequency, and 
A is the sine-wave amplitude. 

* Notice that Slepian's P[T, r(r)] equals one half of our Q(T). 
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As derived in the next subsection, the result obtained is 

Q(T) ^ q? 1 + g< 2) 



where 



with 



q' 2 2> - 2 f (1 - u)ff(i*D 

•'0 



(/(/ 
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(25) 
(26) 



q /.lin-' p(uT) ( 

H(uT) =^J exp<- 



A 2 1 — sin 6 cos 27t&m \ 
2 " cos 2 5 J 

'A 2 si n fl — cos 2ir bu 
cos 2 



dfl 



b = jT 



(27) 



and /„(x) = modified Bessel function of the first kind, zero order. 

The expression given below for gj" is approximate, except for T = k/f, 
k — 0, 1, 2, • • • where it is exact and consequently the function is most 
accurate in a neighborhood of these points. In addition, the accuracy 
of the approximation improves as A increases. For small A, where the 
approximation is least accurate, (26) dominates q^ and so very little 
error results in the upper bound, (25) for all values of A . The expression 
for q ( 2 l) is given by 

• -.(.-» 

/ (b -n - v)S(v) dv, n ^ b <n + \ 

Jo 

- I (6 - n - \ + v)S(v) dv 
Jn+i-b 

+ f\b -n - v)S(v) dv, n + \^b <n + \ 

Jo 

pb- n -i 

- {b -n - \ - v)S(v) dv (28) 

•'0 

+ 2 f (i - v)S(v) dv, n + |^6<n + f 

Jo 

f (b - n - 1 + v)S(v) dv 

Jn+l-b 

- J (b -n - 1 + v)S(v) dv, n + l^b <n + l 

J I) 



(1) ~ — s/ 

<?2 = t 2 A 
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where 



n = 0, 1, 2, 

S( V ) = (l - 49) - - - G(v) (29) 



-*-/■ 



T 



and 



G(v) = f~ " e-* M c " ax dx - f" e A ' (r) cnax dx (30) 

A 2 

While (28) appears to be a formidable equation, it turns out to be 
easily computed, partly because S(v) does not depend on b and hence 
needs only to be computed once for each value of A. The expression 
(26) , on the other hand, turns out to be time-consuming to compute, 
particularly for large values of b. 

4.1 Derivation of Upper Bound, Given by (25) 

The expression for q>, (14), with x(t) specified by (24) is 

92 = h / / R{ti ' t2) dti dt * (32) 

with R(t lf t 2 ) given by (15). Notice that the expectation in this case 
ranges over the three random variables, w{ti), w(t 2 ) and <p. For con- 
venience, define 

Wx ,t 2 )=~l K'. , k) d<p (33) 

Zir J_ T 

where 

r(ti , t 2 ) = E {sgn [to, + a,] sgn [w 2 + a 2 ]\ (34) 

and 

Wi = w{ti) 

a, = A cos (2vfti + <p) i = 1, 2. 

The latter expectation is with respect to w x and w 2 only. Writing out 
(34) in terms of the definition of E{ • } results in 

r(«i , t 2 ) = t; 2TvTl / / sgn [w, + a,] sgn [w 2 + a 2 ] 

Jr\l — P \T)\ J- x J-„ 

/ w\ + wi - 2 P ( T )w 1 w 2 \ 
■ exp < - 2 r t _ 2/ n j dw y dw 2 (3o) 

where t = t> — ti . 
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Applying Price's theorem, 10 to this equation results in 

3r (*' > *») 2 I g\ + g 2 2 - 2p( T)a,q 2 \ 

dp(r) " tt[1 - p'(r)]' 6XP I 2[1 - P \r)] J" (36) 

Integrating (36) and applying the appropriate boundary condition 
yields 

>•(/,,/,) ='^' / 2 )lp<^o + / o P(T> ^^rfp 

or, 
/(/, , t 2 ) = 4 erf (a,) erf (a 2 ) 

2 f T} 1 / af + a 2 — 2aa,a 2 \ , ,„_. 

+ ; J. a"=T5» p i — 2[i - ar I *" (3/) 

where 

As a result of the natural separation of (37) into the sum of two 
quantities, define 

ff 2 = «{" + rf" (38) 

where, by substituting (33) and (37) into (32), the terms in (38) 
may be defined as 

2 r T r T I" 
ql = ~mS I / / erf (a,) erf (a*) cdp dU dt 2 (39) 

„<*> _ _i [' r r r (r) _i_ 

92 ~ir 2 r 2 i ./„ LJ (1 "-«*)* 

I a J + a 2 — 2aa,a 2 \ , , , 

* exp i — n^? — / da d,f dti du - ■ (4 °) 

The detailed steps of simplifying (39) and (40) are relegated to 
Appendices A and B, respectively. In Appendix A we discuss the na- 
ture of the approximation made in arriving at (28) . 

Before applying the results just obtained to specific situations, a 
power series representation for (32) will be given. The series may be 
derived from (39) and (40) by expanding the integrands of these 
functions in their respective Taylor series, evaluating the resulting 
terms and adding the expansions for (39) and (40) together. This 
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procedure results in 

q 2 =-\i (1 - it) sin" 1 P (Tu) du 

If Wo 



£ f a-u)\co S 2r1Tu- P (Tu)] du A 

4.2 Numerical Results and Comparisons 

We evaluated the inequality (25) with the aid of a digital computer. 
For the first case considered, p(r) = e~ M , / = 2 Hz, T ranged between 
and 3.5 seconds and A, the sine-wave amplitude, was either 1 or 10. 
The results of these computations are plotted in Fig. 1. 

Let us first discuss the A = 10 case. The quantity ql u {T) exhibits 
a damped oscillatory behavior much like a plot of [sin (irfT)/(irfT)] 2 
while the quantity q?\T) decays toward zero quite smoothly. For 




2fT 

Fig. 1 — Upper bound of Q(T) for zero crossings of sine wave plus gaussian 
noise. 
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large A , the former term dominates except at its zeros located at T = k/f, 
k = 1, 2, • • • . The values of q 2 {T) are shown as dotted lines on Fig. 1. 
Since Q(T) is a monotonically decreasing function, its upper bound 
can be constructed by drawing horizontal lines between the local 
minima of q 2 (T) and the first intersection of this line with q 2 {T) to 
the right of the minimum. This accounts for the step-like curve drawn 
in Fig. 1 representing the upper bound for Q(T) with A = 10. 

The curve representing .4 = 1, while not exhibiting as fast a decay 
as the curve for A = 10, shows some interesting features. As con- 
trasted to the last case, the q 2 2 \T) term dominates the q 2 x) (T) term 
and consequently much of the oscillatory behavior noted earlier has 
disappeared. 

Another interesting observation can be made when the A = 1 curve 
is compared with the A = curve (gaussian noise alone) for which 
Q(T) is known to equal (2/tt) sin" 1 (e~ T ) when p (t) = e" |rl . (See Ref- 
erence 5.) Notice that for < T < 0.25 the A = 1 curve lies above 
the A = curve while for 0.25 < T < 0.75 the reverse is true. 

This result can be explained by recalling that T = 0.25 represents 
one-half the period of cos (4irt). For intervals shorter than this, the 
sine wave is not likely to cross zero and the effect is to cause fewer 
zero crossings than would be obtained if the sine wave were absent. 
Conversely, for time intervals longer than one-half the period (T = 0.25 
in this case), the sine wave is sure to cross zero and therefore tend to 
increase the number of zero crossings over the noise-alone case. 

As a result of this observation, it seems reasonable to conjecture 
that for T greater than one half the sine-wave period Q{T), for noise 
alone, also forms an upper bound to Q(T) for the sum of a sine wave 
plus noise. 

Additional calculations were made for comparison with Cobb's 
previously reported approximate results. 1 The quantity that Cobb 
derived is an approximate expression for the probability distribution 
function of zero-crossing intervals, denoted by P (T). Rice gives the 
relationship between Q(T) and P {T) in Reference 4 as 

Q(T) - 1 - 2vT + 2v I f P (t) dt dx (41) 

where v = expected number of zero crossings of (24) . 
As observed in Fig. 1 and 2 of Reference 1, 
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for the large sine-wave amplitudes where Cobb's approximation is 
valid. Thus, 



(42) 



(43) 



»2fT ~v 

Q(T) 9* 1 - 2/T + / / PM <1s dy. 

Jo Jo 

Cobb shows in equation 52 of Reference 1 that 

p ° (s) - ^ exp r i? 5 ~J 

where 

[2(1 +PJ] 4 

P. = P (2fT). 

The approximation (43) is only valid for a « 1. 
Substituting (43) into (42), we obtain 

Q(T) ex (i - 2/T)[o.5 + erf ( L ^ / -) J 

^{-[-|( L ^ a: )']--(-^ » 

As in Reference 1, set 

p(r) = ^ (4.")) 

A = 3 (40) 

2tt/ = 0.875 rad/sec. (47) 

Figure 2 compares the approximate solution based on Cobb's re- 
sults (44), and our upper bound (25). For 2fT < 1, the approximate 
solution is somewhat smaller than the upper bound. For 2/T > 1, the 
approximation becomes negative and therefore of little interest while 
the upper bound gradually approaches zero as T increases. 

V. EXTENSIONS TO OTHER CASES 

The specific applications discussed should not be considered ex- 
haustive. For example, the case where x{t) is the sum of a sine wave 
plus gaussian noise could easily be extended to x(t) being the sum 
of a square wave plus gaussian noise. Although the specific formulae 
may be more complex, the general result (equations 12 and 14) is 
still applicable for x(t) nonstationary or nongaussian. 
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Fig. 2 — Comparison of upper bound and approximate solution of Q(T) for 
zero crossings of sine wave plus gaussian noise. 

In addition, the derivation of the general result can be modified 
slightly to obtain a useful upper bound on the conditional probability 
that x(t) does not cross the zero axis for an interval of length T, 
given that x(t) = at the start of the interval. Slepian has inten- 
sively investigated this latter probability. See Reference 5 for his 
discussion. 
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APPENDIX A 

Derivation of 52 " 

We seek a simpler expression for the term 



?2 U = -y* / / / erf (a,) erf (a 2 ) d<p dt t 



dt 2 



(48) 
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which was first encountered in (39) and where 

a. = A cos (2wft + <p), i = 1, 2. (49) 

Substitution of the definition of erf (x) results in 

gi" -nsmf f Pik , k) dt, dt 2 (50) 

7T 1 Jq Jo 

where 

Pit, , Q = f f [' exp [-\{x 2 + i/ 2 )] (to dy tf*>. (51) 

The first step is to notice the following three easily established 
properties of (51) : 

Pit, , t 2 ) = Pit, - Q = Pit, - (52) 

P(r) = p(t + ^) , n = ±l,±2,..- (53) 

4 +i)- -*&-*)■ (54) 

As a result of (52) , (50) may be written as 

tf ' =~L C iT- r)Pir) dr. (55) 

It is a simple matter to demonstrate that, for any function, H (t) , 
satisfying the requirements of (52) through (54) , 

/ (T ~ r)H(r) dr = 0. (56) 

h/t 

We next introduce an approximation to (51) that preserves properties 
(52) through (54). It is important to preserve these properties because, 
as a result of (56), if they are satisfied, q™ = for T = k/f, fc - 1, 2, • • • ; 
consequently, an approximation satisfying (52) through (54) will be 
accurate in a vicinity of these values of T. In addition, the three prop- 
erties permit fast computation of (50). 

The approximation chosen is given by 



p f" exp [-\ix 2 + y 2 )] dxdy^ jT JT e~ r ' /2 r dr dc 



(57) 
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for sgn [a,] = sgn [a 2 ], and 

I I exp[-$(x+y)]dxdy^-j' j \~ r ' /2 r dr da (58) 

for sgn [ai] = -sgn [a 2 ]. 

Essentially the approximation results in deforming the region of 
integration, as shown in Fig. 3. From this figure, it may be noticed that 
(57) is in reality an upper bound while (58) is a lower bound. Of course, 
it is easy to conceive of functions that give an upper bound to (58) 
and thus result in an upper bound for q?\ However, this results in a 
loss of properties (52) through (54). 

Evaluating the integrals appearing in (57) and (58) and then sub- 
stituting into (51) yields 



Hu ,u) - 1 f' p(fP, u , fe)[i - e -* ( -' +o -'>] d* 



where 
and 



p(<p, h , Q = sgn [a, aj 



(59) 



(60) 



P{U , t 2 ) &£ P(t x , t 2 ). 

Proving that (59) possesses the properties (52) through (54) only 
requires the use of elementary integration theory and will therefore 
be omitted. As a result of these properties, (59) may be written as 



where co = 2rf 

p(6, t) = sgn [cos 6 cos (cor + 0)]; 

y 



do 



(61) 
(62) 



REGION OF INTEGRATION 
FOR APPROXIMATION 



\ REQUIRED REGION 

^Y" OF INTEGRATION 

\ 

\ , 

^ — radius = i/a, 2 +a 2 2 

> 



Fig. 3 — Deformation of region of integration for approximation in Appendix A. 
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furthermore, P(t) need only be evaluated for the range < r < 1/(4/). 
Values beyond this range are related to values within the range by 
(52) through (54). To evaluate (61), an explicit expression for (62) is 
required. After studying this latter equation one finds that for 

< r < 1/2/ 
1 for I < d^ | - cor 



p(d, t) = , 



1 for | - cot < % | 



1 for - \ < 6 ^ | - cor 



-1 for - I - cor < 9 ^ -| 



(63) 



with similar expressions for t falling in the ranges 

k = 1,2, • 



2/ = - 2/ 



However, only the expression given by (63) is needed to evaluate 
(61) in the required range. 
Substitution of (63) into (61) results in 

U(»/2)-o,r «CW2)- W r 

F(r, d)dd+ F(t, 6) dd 

■t/2 J */2 



- / F(t, 6) dd - / F(t, 6) dd 



(64) 



where 



F(t, d) = 1 - e 



-(AV2)[oos'J+ cos' <wr + < 



With the help of some fundamental trigonometric identities it is easy 
to show that 



F(r, d) = 1 - e- At "e- KW 0O ' (wT+29) 



(65) 



where 



K(t) = -jr cos cor. 



(66) 
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Substituting (65) into (64) and performing obvious simplifications 
results in 

P(t) = I {2tt - 4cor - 2e- A ' /2 G(r)\ (67) 

where 

G( T ) = f "'«-«« — &- rV w — lb. (68) 

Jo Jo 

For the time being assume ra ^ jT < n + J where n = 0, 1, 2, • • • . 

Substituting (67) into (55) gives 



ffj" S p jT (T - r)fl(r)di 



where 

5(r) = 1 - 4/r - ? — 0(r). (69) 

7T 

Using the result (56) , the latter equation equals 

tf ' ^~ f (T- r)S(r) dr. 
Setting t = t — n/f, 

~ 5* I \ T -J- ( >«> * 

as a result of property (53). Now substitute t = v/f to obtain 

q? = |s / (b-n- v)S(v) dv for n £ b < n + J (70) 

where 6 = /7 1 . 

This equation is the same as the first part of the final result stated 
in (28). The equations defined in (69), (68), and (66) are the same 
as (29), and (30), and (31), respectively, except for a convenient 
scale change. The rest of the results stated in (28) , for various 
ranges of T, are derived in a similar manner as (70) was, using rela- 
tions (52) , (53) , or (54) , as required. Because only straightforward 
operations are used to obtain these results, they will not be derived. 
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APPENDIX B 

Derivation of q {2) 

The first step in simplifying the expression for q 2 2) , as defined in (40), 
is to interchange the order of integration of the two innermost integrals 
to yield 

ft™ = ~Ls f f f () n X *\i B( > a > l * ' ^ da dtl dU - (71) 

IT 1 Jo Jo Jo \1 a ) 

where 



B(a, ti ,t 2 )=\ exp 



a 2 + a l — 2aa l a 2 \ , n9 . 

2(1 - « a ) J *" (?2) 

Using the definitions of a x and a 2 , (34), and some obvious trigono- 
metric identities, it is easy to show that 

a\ + a\ — 2aa l a 2 

= A 2 {cos [w(/, + t 2 ) + 2<p][a - cos («r)] + [1 - a cos («t)]| 

where we have set w = 27r/ for convenience. 

Substitution of this relationship into (72) results in 

B(a, ty , t 2 ) = exp [-</,(«, t)] 

•/'exp (-/,(«, r) cos [»(*, + k) + 2^]} dip (73) 



where 



_ . . A 2 f l — « cos (cot) 



t , \ A? [ a - COS (ot) 
J 2 (a, r) = y\ 1 _ a » 



(74) 
(75) 



Setting = ta(k + fc) + 2*> in (73) and using the periodic properties 
of the integrand, yields 



B(a, r) = 2 exp [-,/,(<*, t)] / exp [-«/ 2 (a, t) cos 0] 

Jo 



dd. 



This integral is recognized as an expression for the modified Bessel 
function of the first kind (see Reference 11, page 181, Equation 4). 
And so 



B(a, t) - 2irI [J 2 (a, r)] exp [-Ji(«, r)] 



(76) 
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where I (x) is the modified Bessel function of the first kind, zero 
order. 

Since (76) is only a function of r, one may define 

ff to = - f ( ° 7= — 1 *(«» <**■ ( 7? ) 

7T J VI- a 

Furthermore, it is easy to show that H(t) = H(—r). Consequently, 
(71) can be written as 



qi 2) = 2 f (1 - u)H(uT) du. 

Jo 



By setting t = uT in (77) and by making the change of variable a 
= sin 0, (26) and (27), which are the desired results, follow. 

REFERENCES 

1. Cobb, S. M., "The Distribution of Intervals Between Zero Crossings of Sine 

Wave Plus Random Noise and Allied Topics," IEEE Trans. Inform. 
Theory, 1T-11, No. 2 (April 1965). pp. 220-231. 

2. Newell, G. F., and Rosenblatt, M., "Zero-Crossing Probabilities for Gaus- 

sian Stationary Processes," Ann. Math. Stat., 33, No. 4 (December 1962), 
pp. 1306-1313. 

3. Rice, S. 0., "Mathematical Analysis of Random Noise," B.S.TJ., S3, No. 3 

(July 1944), pp. 282-332; and 2.',, No. 1 (January 1945), pp 46-156. See 
especially Sections 3.3 and 3.4. 

4. Rice, S. O., "Distribution of the Duration of Fades in Radio Transmission," 

B.S.T.J., 37, No. 3 (May 1958), pp. 581-636. 

5. Slepian, D., "The One-Sided Barrier Problem for Gaussian Noise," B.S.T.J., 

41, No. 2 (March 1962), pp. 462-501. 

6. Wong, E., "Some Results Concerning the Zero-Crossings of Gaussian Noise," 

SIAM J. Appl. Math., 14, No. 6 (November 1966), pp. 1246 T 1254. 

7. McFadden, J. A., "Urn Models of Correlation and a Comparison with the 

Multivariate Normal Integral," Ann. Math. Stat., 26, No. 3 (September 
1955), pp. 478-489. 

8. Strakhov, N. A., "A Representation Theory for a Class of Non-Gaussian 

Distributions with Applications to Digital Data Transmission," New 
York: Doctoral Dissertation, New York University, October 1967. 

9. Lawson, J. L., and Uhlenbeck, G. E., Threshold Signals, New York: Mc- 

Graw-Hill Book Co., Inc., 1950. 

10. Price, R., "A Useful Theorem for Nonlinear Devices Having Gaussian In- 

puts," IRE Trans. Inform. Theory, 1T-4, No. 2 (June 1958), pp. 69-72. 

11. Watson, G. N., A Treatise on the Theory of Bessel Functions, 2nd ed., 

Cambridge, England : University Press, 1962. 



